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Graphene nanodisk is a graphene derivative with a closed edge. The trigonal zigzag nanodisk with size N has 
iV-fold degenerated zero-energy states. A nanodisk can be interpletted as a quantum dot with an internal degree 
of freedom. The grand state of nanodisk has been argued to be a quasi-ferromagnet, which is a ferromagnetic- 
like states with a finite but very long life time. We investigate the spin-filter effects in the system made of 
nanodisks and leads based on the master equation. The finite- size effect on spin filter is intriguing due to 
a reaction from the polarization of incoming current to a quasi-ferromagnet. Analyzing the relaxation process 
with the use of the Landau-Lifshitz-Gilbert equation, we explore the response to four types of incoming currents, 
namely, unpolarized current, perfectly polarized current, partially polarized current and pulse polarized current. 
We propose some applications for spintronics, such as spin memory, spin amplifier, spin valve, spin-field-effect 
transistor and spin diode. 



I. INTRODUCTION 

The study of spin-dependent transport phenomena has re- 
cently attracted much attention. 1,2 It has opened the way for 
the field of spintronics, 3,4,5,6,7 literally spin electronics, where 
new device functionalities exploit both the charge and spin de- 
grees of freedom. There are various approaches in this sphere. 
For instance, the use of a quantum-dot setup 8,9 has been pro- 
posed, which can be operated either as a spin filter to pro- 
duce spin-polarized currents or as a device to detect and ma- 
nipulate spin states. Spintronics in graphene has also been 
proposed 1 ! 11 ! 12 recently. In the following we address this is- 
sue in a new type of materials, graphene nanodisks. 1 3 1 14 1 15 1 16 1 17 

Graphene nanostructure 18,19,20 such as graphene 
nanoribbono 21 i 22 i 23 i 24 i 25 i 26 i 27 i 28 i 29 i 30 and graphene 

nanodisks 13 i 14 i 15 i 16 i 17 could be basic components of fu- 
ture nanoelectronic and spintronic devices. Graphene 
nanodisks are nanometer- scale disk-like materials which have 
closed edges. They are constructed by connecting several 
benzenes, some of which have already been manufactured 
by soft-landing mass spectrometry ] 27 i 31 i 32 There are varieties 
of nanodisks, among which zigzag trigonal nanodisks with 
size N are prominent because they have iV-fold degenerated 
zero-energy states. We have already shown 14 that spins make 
a ferromagnetic order and that the relaxation time is finite 
but quite large even if the size N is very small. We refer to 
this property as quasi-ferromagnet. Furthermore we have 
argued 17 that a nanodisk behaves as if it were a quantum dot 
with an internal degree of freedom, where the conductance 
exhibits a peculiar series of Coulomb blockade peaks. 

In this paper we make an investigation of the spin current 
in the zero-energy sector of the trigonal zigzag nanodisk. We 
first analyze how the spin of a nanodisk filters the spin of 
the current by assuming that the nanodisk is a rigid ferro- 
magnet. The fact that the direct and exchange Coulomb in- 
teractions are of the same of magnitude plays an important 
role. However the nanodisk is not a rigid ferromagnet but a 
quasi-ferromagnet. Hence an intriguing reaction phenomenon 
occurs: Namely, the spin of the nanodisk can be controlled 
by the spin of the current. Using these properties we propose 
some applications for spintronic devices, such as spin mem- 



ory, spin amplifier, spin valve, spin-field-effect transistor and 
spin diode. 

This paper is organized as follows. In Sec|IIl we summarize 
the basic notion of the nanodisk. The low-energy physics is 
described by electrons in the TV-fold degenerated zero-energy 
sector, which form a quasi-ferromagnet due to a large ex- 
change interaction. In Sec|IIIl we investigate the spin-filter 
effects based on the master equation. In Sec|IVl we analyze 
the reaction to the spin of the nanodisk from the spin of elec- 
trons in the current. We discuss the relaxation process using 
the Landau-Lifshitz-Gilbert equation. In Sec|Vl we propose 
various spintronic devices and explore their properties. 



II. NANODISK QUASI-FERROMAGNETS 
A. Zero-Energy Sector 

Graphene nanodisks are graphene derivatives which have 
closed edges. The Hamiltonian is defined by 



H 



(ho) 



(2.1) 



where e\ is the site energy, tij is the transfer energy, and c\ 
is the creation operator of the tt electron at the site i. The 
summation is taken over all nearest neighboring sites 
Owing to their homogeneous geometrical configuration, we 
may take constant values for these energies, Si = and Uj = 
t « 2.70eV. 

We have searched for zero-energy states or equivalently 
metallic states in graphene nanodisks with various sizes and 
shapes. 14 The emergence of zero-energy states is surprisingly 
rare. Among typical nanodisks, only trigonal zigzag nan- 
odisks have degenerate zero-energy states and show metallic 
ferromagnetism, where the degeneracy can be controlled arbi- 
trarily by designing the size. 

Trigonal zigzag nanodisks are specified by the size param- 
eter N as in FigQIa). The size-TV nanodisk has TV-fold degen- 
erated zero-energy states, 14 where the gap energy is as large as 
a few eV. Hence it is a good approximation to investigate the 
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electron-electron interaction physics only in the zero-energy 
sector, by projecting the system to the subspace made of those 
zero-energy states. The zero-energy sector consists of N or- 
thonormal states \f a ), ol = 1, 2, • • • , TV, together with the 
SU(N) symmetry. We can expand the wave function of the 
state \f a ), 

fa(x)=Y^LJ?<Pi(x), (2.2) 

i 

in terms of the Wannier function (fi(x) localized at the site i. 
The amplitude wf is calculable by diagonalizing the Hamilto- 
nian (12.11) . 



(a) nanodisk (b) nanodisk-lead system 




FIG. 1: (a) Trigonal zigzag nanodisks. The size parameter N is 
defined in this way. The number of carbon atoms is given by 
Nc = N 2 + 6N + 6. (b) The nanodisk-lead system. The nanodisk 
with N = 7 is connected to the right and left leads by tunneling 
coupling £r and £l- The chemical potential is /xr or /xl at each lead. 



B. Quasi-Ferromagnets 

We include the Coulomb interaction between electrons in 
the zero-energy sector. 17 It is straightforward to rewrite the 
Coulomb Hamiltonian as 

Hf = Y, U «pn{a)n{(3) 

a>p 

~ \ J2 J^S(a) • S((3) + n (a) n (/?)], (2.3) 

a>p 

where U a p and J a p are the Coulomb energy and the exchange 
energy between electrons in the states \f a ) and \fp). Here, 
n (a) is the number operator and S(a) is the spin operator, 

n (a) = dl(a)d a (a), (2.4a) 

S(a) = ^4(«) r ^d ff /(a), (2.4b) 

with d<j{a) the annihilation operator of electron with spin 
a = j, | in the state \f a )'- t is the Pauli matrix. 

The remarkable feature is that there exists a large overlap 
between the wave functions f a (x) and fp(x), a ^ (3, since 
the state \f a ) is an ensemble of sites as in (12.21) and identical 
sites are included in \f a ) and \fp). Consequently, the domi- 
nant contributions come from the on-site Coulomb terms not 
only for the Coulomb energy but also for the exchange energy: 
Indeed, it follows that U a p = J a p in the on-site approxima- 
tion. We thus obtain 

U a p~J a p~U^(LJ?LJ?) 2 , (2.5) 



where 

U = J2 [d 3 xd 3 ytf(x)^(x)V(x-y)^(y)^(y) (2.6) 

with the Coulomb potential V(x — y). The Coulomb energy 
U is of the order of leV because the lattice spacing of the 
carbon atoms is ~ 1A in graphene. 

Since the exchange energy J a p is as large as the Coulomb 
energy U a p, the spin stiffness J a p is quite large. Furthermore, 
we have checked 17 numerically that all J a p are of the same 
order of magnitude for any pair of a and (3, implying that the 
SU(N) symmetry is broken but not so strongly in the Hamil- 
tonian (12.31) . It is a good approximation to start with the exact 
SU(7V) symmetry. Then, the zero-energy sector is described 
by the SU(N) Heisenberg-Hubbard model, 

H D = -J^S(a)-S(f3) 

a^p 




with J ~ U . We rewrite as 




(2.8) 



in terms of the total spin, 

S toi = Y,S(a), (2.9) 

a. 

and the total electron number, 

ritot = (2.10) 

a 

At the half-filling, the eigenstate of the Hamiltonian 
is labeled as \s) = \n tot ,l,m), where I is the total spin- 
angular momentum and m is its z-component. Based on this 
Hamiltonian we have shown 14 that the relaxation time of the 
ferromagnetic-like spin polarization is quite large even if the 
size of trigonal zigzag nanodisks is quite small. Such a nan- 
odisk may be called a quasi-ferromagnet. We refer to the total 
spin S t ot of a nanodisk as the nanodisk spin. 

C. Nanodisk-Lead System 

We consider a system made of a nanodisk connected by two 
metallic leads [FigQJb)]. The model Hamiltonian is given by 

H = H L + H TL + H TR + H D , (2.11) 
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where H D is the Hamiltonian (12.7b of a nanodisk, and H L is 
the lead Hamiltonian 

H h = (k) (c£ct + c^cL) • (2.12) 

ka 

The Hamiltonian i^L describes a noninteracting electron gas 
in the leads with e (k) = h 2 k 2 /2m, while ifxL and Hjr are 
the transfer Hamiltonian between the left (L) and right (R) 
leads and the nanodisk, respectively, 

#tl =*l E E ( C fel^(«) + 4(«)cL) , (2.13a) 

ka ol 

H TR =t R E E ( C *K(") + 4(«)cL) , (2.13b) 

/ecr a 

with t x the tunneling coupling constant: We have assumed 
that the spin does not flip in the tunneling process. 

The nanodisk-lead system looks similar to that of the TV-dot 
system. 33 However, there exists a crucial difference. On one 
hand, in the ordinary TV-dot system, an electron hops from one 
dot to another dot. On the other hand, in our nanodisk system, 
the index a of the Hamiltonian runs over the TV-fold degen- 
erate states and not over the sites. According to the Hamilto- 
nian (12.131) . an electron does not hop from one state to another 
state. Hence, it is more appropriate to regard our nanodisk as 
a one-dot system with an internal degree of freedom. 

III. SPIN FILTER 

The aim of this paper is to investigate a nanodisk as a spin 
filter 38 based on the Hamiltonian (12.111) . The setup we con- 
sider is a lead-nanodisk-lead system [FigQIb)], where an elec- 
tron makes a tunnelling from the left lead to the nanodisk and 
then to the right lead. We investigate how the electron spin is 
affected by the nanodisk spin during the transport process. 

The dynamics of the nanodisk system is described by the 
master equation, 

^ = E [W (s, s') p (s>) - W (s>, s) p (s)j , (3.1) 

s' 

where p (s) represents the probability to find the system in the 
state \s) = |n to t, m), and W (V, s) is the transition rate be- 
tween the states \s) and \s'). The master equation describes 
a stochastic evolution in the space spanned by the states \s). 
Jumps between different states are assumed to be Markovian. 
The stationary solution is given by the detailed balance condi- 
tion dp/dt = 0. 

Since the Hamiltonian H D only depends on S? oV the proba- 
bility p(s) has no dependence on the spin component m. Then 
it is convenient to do a sum over m, 

P(s) = ^2p(s) (3.2) 

m 

with \s) = |n to t, 0- The master equation is rewritten as 

^P=£[^ & P&)-W (s', s) -p (!)] . (3.3) 



In the following we denote p (s) as p (s) for simplicity. 

In this paper we consider the small coupling limit, t x <C 
J, x =L, R, where the dominant process is the sequential 
tunneling: It is of the order of |t x | 2 , while the cotunneling 
process of the order of \t x \ . 

In the sequential-tunneling regime the transition rate, W = 
Wl + Wr, is obtained by the Fermi's golden rule, 8 ^*^ 

W x (V, s) =T x f x (A s / S ) £ n /, n+ i 

+ r x [i-/ x (A fl / fl )](y n /, n _i. (3.4) 

Let us explain notations. First, 

f x (e) = {1 + exp [( £ - p x ) /fcBT]}" 1 (3.5) 

is the Fermi function at temperature T, with p x the chemical 
potential at the right or left lead, x =R,L- We set 

Ml = /i + A/i, p R = p- Ap. (3.6) 

Second, A s / S is the energy difference between the two states 
\s'} and \s), 

A s , s =E 8 > -E s . (3.7) 

Third, 

T x = 2w\Al ssf \ 2 (3.8) 

is the tunneling rate through the right or left lead, where v — 
S (ep — £k) is tne density of states at the Fermi level £p in 
the leads, which is a constant, and 

A x , ss ,=t x Y,(s'\d aa \s). (3.9) 

a 

It follows that the tunneling rate T x is a constant for the states 
\s) and \s f ) connected by a Markov step, and zero otherwise. 
The current through the nanodisk can be written as 

la = ^' P ( S ) - ^ R S ') P ^ ' < 3 - 10 ) 

s,s' 

where J2 S s' runs over me states I s ) an( i I s ') sucn mat 
(s'\ d aa \s) 7^ 0. Thus the current I a depend on the spin con- 
figuration of electrons in the nanodisk. 

We have argued that the nanodisk is a quasi-ferromagnet. 
For the sake of simplicity, we start with the simplification 
that the ground state is a ferromagnet with polarized up- spins. 
When the transition interaction is small enough, a single elec- 
tron tunnels at once, namely, the electron number in the nan- 
odisk increases or decreases by one, n[ ot = n tot ± 1. First 
of all, we notice that, when n[ ot = n to t + 1, the tunneling of 
spin-up electrons is blocked by the Pauli's exclusion principle. 
The Markov chain of spin configurations contain only a finite 
number of states. For instance it is 9 in the case of TV = 5 
nanodisk. We number them as indicated in Fig|2j 

We calculate the energy E s of various states \s) = |n to t, I) 
in this chain. According to the above numbering convention 
(Fig|2|), they are 

E\ =E NjN / 2 , E 2 = En-i,(n-i)/2, ^3 = En, N/2-1, 

E4 =E N _ 1 ^ N _ 3 y 2 : E$ = En+1,(N-1)/2, 

(3.11) 
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FIG. 2: Markov chain of spin configurations in the N — 5 nanodisk. 
An arrow (|) indicates each Markov step between two states with 
different spin configurations. The region surrounded by red solid 
lines shows the most dominant process for spin-polarized current. 
In this approximation there is no down- spin polarized current. The 
region surrounded by green dotted lines yields the next dominant 
contribution, which allows a down-spin polarized current as well. 



Based on the effective Hamiltonian J2.8I) they are calculated 
as follows, 
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N 



-N « UN, 



(3.12a) 



(N-l) 2 + ^±A(N-l)*U(N-l), 



N 



-N « 2UN, 



(3.12b) 
(3.12c) 



(n - iy + 



2 /7 + 3J 



(JV-1)«217(JV-1), 
(3.12d) 



N 2 + - (U + J) TV - J * 3/77V - 17, 



and so on. 

The order of the energies is 



E 2 < E x < E 4 < E 3 < E b < 



(3.12e) 



(3.13) 



The probability to find the state \s) contains the Boltzman fac- 
tor exp(— E s /k^T). Hence, the most dominant Markov chain 
consists only of the states |1) and |2), and it is denoted as 
(1 <-> 2). This process yields the spin-up current. We consider 
the second dominant process. Though E4 < E3, the state |4) 
can only be reached via the state |3), as illustrated in Figj2] 
Namely, the second dominant chain is (1 <-> 2 <-» 3 <-> 4), 
which allows the spin-down current as well as the spin-up cur- 
rent. 

First we analyze the most dominant process, which involves 
only the ground state and the first excited state (Figj2j). The 
currents are given by 

I 1 = I -^[f L (E l -E 2 )-U(E l -E 2 )} 



sinh(A/i/fc B T) 



2 cosh(A A t/fc B T) + cosh((C/ - n)/kaT) ' (3 ' 14a) 
h =0, (3.14b) 



where 



er L r R 



(3.15) 



r L + r R 

in terms of the tunneling rate T x given by (13.8b . 

The above expression can be further simplified in the fol- 
lowing two limits. In the high-temperature limit (/t, A/i <C 
k^T) the sequential tunneling current takes a simple form 



4A„ lim 9f ^~ 
2 r a m ^o dAfi 



E 2 ) 



T A I Jj U-2 
=in — COSh 

8/c B T 



*12 



/i 



2/c B T 



(3.16) 



where A s / S is defined in (13.7b . On the other hand, in the low- 
temperature limit (k^T <C /i, A/i) the current is given by 

T 2 
with the step function 



[0 (/i L - A 12 ) - (/^ - A 12 ) 



(3.17) 



>(*) = 



x > 
x < 



(3.18) 



The current 7| flows in the triangle domain (/i^ > A12 > /ii?) 
in the (A/i, /i) plain. 

In order to search for the down- spin polarized current 7j, 
it is necessary to analyze the second dominant process. It in- 
cludes four states (1 <-> 2 <-> 3 <-> 4), as shown in the Fig 12 
In this approximation, the transition matrix becomes a tridi- 
agonal matrix, and we can solve the master equations exactly. 
The nonequilibrium density of states is calculated as 



Trp s 



Ps 



n 



W(i + l,i) 

\ w(i,i + iy 



(3.19) 



where the index s = 1, 2, 3, 4 stands for the four states with 
the energy E s . We can explicitly calculate p s as 



-0E, 



n 



e /3B i+ i cos h A/i 



\ eP Ei cosh A/x + e^+i ' 



(3.20) 



with E s = E s — /intot- The current /,j is calculable with the 
use of (13.101) and (13.41) . whose result we show in the Fig|3] In 
particular we have 



x [9 (ml 



A 12 ) 
A 32 ) 



A 12 ) 
A32) 



(3.21) 



in the zero -temperature limit. The current 1^ flows in the tri- 
angle domain (/i L > A i2 > Hr and \±l> A 32 > \ir) in the 
(A/i, /i) plain. 



IV. FINITE-SIZE EFFECTS: REACTION TO 
QUASI-FERROMAGNET 

In the conventional spin filter the ferromagnet is very rigid. 
On the other hand, the life time is finite in the case of a nan- 
odisk, though it is quite long in spite of its small size. We 
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FIG. 3: Spin-dependent currents against n/k^T and A/x/ZcbT in the 
N — 10 nanodisk. We have set U '/kaT — 1. The horizontal x-axis 
is the chemical potential fi/kaT, and the horizontal y-axis is the bias 
voltage A/x//cbT. (a) The vertical z-axis is the current, (b) More 
currents flow in brighter regions. (1) The current I\ induced by the 
process 1 2, which produces the up-spin polarized current. (2) 
The current h induced by the process 3^4, which produces the 
up-spin polarized current. (3) The sum of the above two processes; 
ij = 1 1 + h. (4) The current I a induced by the process 2^3, 
which produces the down- spin polarized current; i| = I a. (5) The 
difference A I between the up-spin and down- spin current, A I — 
J T - I i = h + h - h. 



expect a reaction to the nanodisk spin from the spin-polarized 
current, provided the nanodisk spin is not controlled exter- 
nally, say, by magnetic field. This effect is very interesting, 
since it is an intrinsic nature to quasi-ferromagnets. 

We inject an electron into a nanodisk. Let A^ or A^ be 
the energy increase when the spin of an injected electron is 
parallel or anti-parallel to the nanodisk spin. According to the 
Hamiltonian (12.71) they are 



Att 



A 2i =-U, 
A 3 i =UN. 



(4.1) 



When the direction between the nanodisk spin and the electron 
spin is 9, the energy increase is given by 

9 6 

(Hg f ) = A TT cos -(c{c T ) + A n sin -<c{q), (4.2) 

where {c^c a ) represents the probability of finding the injected 
electron to have spin a. 

Based on the spin-rotational symmetry we may write the 



effective Hamiltonian as 

h d = ~ J sdS ■ c ] ac = -Mn ■ c ] ac, (4.3) 

where 

J sd = MS tot = A n — Aft = U(N + 1), (4.4) 

and ft — Si i/ Sxqx is the normalized spin with /St t — l^totl- 

We introduce the Gilbert damping term phenomenologi- 
cally, 



H a = an , 



(4.5) 



where a is a dimensionless constant (a « 0.01). Using the 
variational method for the Hamiltonian H^ f + H a , we obtain 
the Landau-Lifshitz-Gilbert equation, 



On On 

— jB Q ff xn-anx — , 
at at 



(4.6) 



where £? e ff is the effective magnetic field produced by the in- 
jected electron spin, 

with the gyromagnetic ratio 7. 

For definiteness we now inject the up-spin electron to the 
nanodisk. The effective magnetic field is 



7 #eff= (0,0, 7 5 eff ) 



M 

hStot 



(cW) . (4.8) 



We investigate the dynamics of the normalized spin n of the 
nanodisk under this field. In the polar coordinate, setting 

n = (sin cos 0, sin sin 0, cos 0) , (4.9) 

we rewrite the Landau-Lifshitz-Gilbert equation (14.61) as 

= a0sin0, </>sin0 = 7# e ff sin0 - a0. (4.10) 

This is equivalent to 



a^B { 



eff 



■ sm e/, 



7^eff 

1 w 



1 + a 2 

which we can solve explicitly, 

(t) = 2 tan -1 exp [- (t - t ) /r filte r] . 

1 + cr 



(4.11) 



(4.12) 



where to is an integration constant, and Tfii ter is the relaxation 
time given by 



Tfiiter OC N. 



2a7|,Beff| 



(4.13) 



It is proportional to S tot /M and hence proportional to the disk 
size TV because of (14.71) . The initial phase 0o is related to the 
parameter to as 



(0) = 2 tan 1 exp [t /r mteT \ . 



(4.14) 
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FIG. 4: The bold curve is the relaxation process of spin in nan- 
odisk. The horizontal axis is the time t/r and the vertical axis 
is the angle 0. The thin curve is the asymptotic function = 
2exp [— (t — to) /^filter] of the bold curve for t ^> to- The dotted 
line is the asymptotic value for t <C to. 



Thus, to = corresponds to 



'<o) = f. 



(4.15) 



We use the parameter to instead of the initial phase 6(0) for 
simplicity. 

The time scale for the direction of the nanodisk spin to align 
with that of the spin-polarized current is Tfii ter , where the effec- 
tive magnetic field (14.71) is proportional to the injected current 
J in . In other words, we can control the polarization of the 
nanodisk spin by using the spin-polarized current. We note 
that this is possible since the nanodisk is a quasi-ferromagnet. 
Indeed, there exist no effective magnetic field (14.71) in the con- 
ventional ferromagnet. 



V. SPINTRONIC DEVICES AND APPLICATIONS 

We summarize the spin properties of a nanodisk and 
an incoming electric current. First of all, being a quasi- 
ferromagnet, the nanodisk has a definite polarization. With 
respect to the incoming electric current there are three cases. 
(1) The polarized current, where all electrons have a definite 
polarization, rotates the nanodisk spin to that of the incom- 
ing current, as we have shown in SeclIVl (2) The unpolarized 
current, where the polarization of each electron is completely 
random, does not induce any effective magnetic field. Hence 
it is filtered so that the outgoing current is polarized to that of 
the nanodisk. (3) The partially polarized current, where the 
polarization of each electron is at random but the averaged 
polarization has a definite direction, induces a net effective 
magnetic field. Hence it rotates the nanodisk spin to that of 
the incoming current, and then is filtered so that the outgoing 
current is completely polarized to the averaged polarization 
of the incoming current. Furthermore, it is possible to control 
the nanodisk spin externally by applying magnetic field. Then 
the outgoing current has the same polarization as that of the 
nanodisk, irrespective of the type of incoming current. Using 
these properties we propose some applications of graphene 
nanodisks for spintronic devices. 



A. Spin Memory 

The first example is a spin memory. 8 For a good memory 
device three conditions are necessary: (i) It keeps a long life 
time information; (ii) Information stored in the memory can 
be read out without changing the information stored; (iii) It is 
possible to change the information arbitrarily. 

First, since the life time of the nanodisk quasi-ferromagnet 
is very long compared to the size^, 



Tferro CX N 



(5.1) 



we may use the nanodisk spin as an information. Next, we can 
read-out this information by applying a spin-unpolarized cur- 
rent. The outgoing current from a nanodisk is spin-polarized 
to the direction of the nanodisk spin. Thus we can obtain the 
information of the nanodisk spin by observing the outgoing 
current. Finally, the direction of the nanodisk spin can be con- 
trolled by applying a spin-polarized current into the nanodisk. 

Thus, the nanodisk spin satisfies the conditions as a mem- 
ory device. The important point is that the size is of the order 
of nanometer, and it is suitable as a nanodevice. 



B. Spin Amplifier 

The second example is a spin amplifier. We inject a 
partially-polarized- spin current, whose average direction we 
take to be up for definiteness. Thus, If > If > 0. On 
the other hand, the direction of the nanodisk spin is arbitrary. 
Since spins in the nanodisk feel an effective magnetic field 
proportional to If — If, they are forced to align with that of 
the partially-polarized- spin current after making damped pre- 
cession. By using (14.121) the time dependence is given by 



7 T (t) = 4 n cos 
I^t) = If sin - 



2 



1 



^1 + exp [-2t/r fi iter] ' 
1 

y/1 +exp [2t/r fi iter]' 



(5.2a) 
(5.2b) 



where we have set to = 0, which means 6(0) = tt/2: See 
Fig 13 The outgoing current is initially given by 



1 



1 



J T (0) = -j=I? , I, (0) = -^=If. (5.3) 

as a function of the incoming current I™. The time scale is 
given by the relaxation time (14.13b . 

After enough time t — to ^> r, all spins in the nanodisk 
take the up direction and hence the outgoing current 7° ut = 
lim^oo I a (t) is the perfectly up-polarized one, 



jout jin 



(5.4) 



Consequently, the small difference If — If is amplified to 
the large current If. The amplification ratio is given by 
If /(If — If), which can be very large. This effect is very 
important because the signal of spin will easily suffer from 
damping by disturbing noise in leads. By amplifying the sig- 
nal we can make circuits which are strong against noises. 
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1 2 3 4 5 
time (f/r) 

FIG. 5: Relaxation process of spin amplifier. The horizontal axis is 
the time t/r and the vertical axis is the current ij and i|. We have 
set If = 0.6^/2, If = OAy/2. The currents saturate after enough 
time (t > 2t for if, t > 5r for ij), and the amplification ratio is 3 
in this example. 



C. Spin Valve and Spin-Field-Effect Transistor 

The third example is a spin valve, or giant magnetoresis- 
tance effect [Fig|6l . 39 i 4Q i 41 We set up a system composed of 
two nanodisks sequentially connected with leads. We apply 
external magnetic field, and control the spin direction of the 
first nanodisk to be 



|0>=cos-|T> + sin-U>. 
and that of the second nanodisk to be 

|0) = |T). 



(5.5) 



(5.6) 



We inject an unpolarized-spin current to the first nanodisk. 
The spin of the lead between the two nanodisks is polarized 
into the direction of \9). Subsequently the current is filtered to 
the up- spin one by the second nanodisk. The outgoing current 
from the second nanodisk is 



7 T 0Ut = ^T cos-. 

T 1 2 



(5.7) 



Since we can arrange the angle externally, we can control 
the magnitude of the up-polarized current from zero to one. 
In this sense the system act as a spin valve. 

The forth example is a spin-field-effect transistor 36 [FiglT). 
We again set up a system composed of two nanodisks sequen- 
tially connected with leads. We now apply the same external 
magnetic field to both these nanodisks, and fix their spin di- 
rection to be up, 



|o) = lt). 



(5.8) 



As an additional setting, we use a lead between the two nan- 
odisks possessing a strong Rashba-type spin-orbit coupling 37 , 



Hr = ^(p x cr y -p y (T x ). 



(5.9) 



Spins make precession while they pass through the lead. The 
spin-rotation angle is given by 4 



(a) spin-valve transistor 
L 




(b) parallel external fields (spin current flows) 




(c) anti-parallel external fields (spin current is blocked) 



II Kill 



FIG. 6: Illustration of spin valve with unpolarized current incoming 
from left. Spins of each nanodisk are controled by the external mag- 
netic field. The unpolarized current is filtered by the left nanodisk, 
and only electrons whose spin direction is the same as the nanodisk 
spin go through the centeral lead, (a) The spin valve is made of two 
nanodisks with the same size, which are connected with leads, (b) 
We apply the external magnetic field in the same direction. The spin- 
polarized outgoing current flows, (c) We apply the external magnetic 
field in the opposite directions. The outgoing current is blocked. 



where m* is the electron effective mass in the lead and L is 
the length of the lead. We can control AO by changing the 
coupling strength A externally by applying electric field. 42 The 
outgoing current from the second nanodisk is 



rout 

T 



Jf cos 



AO 



(5.11) 



Since we can arrange the angle AO by applying electric field 
and control the magnitude of the up-polarized current from 
zero to one, we expect the system acts as a spin-field-effect 
transistor as in the conventional case. 



D. Spin Diode 

The fifth example is a spin diode [FiglS). We use a system 
similar to the spin-field-effect transistor but with the following 
differences. First, two nanodisks have different sizes. When 
the left nanodisk is larger than the right nanodisk, the relax- 
ation time of the left nanodisk tl(= Tj~j ter ) is larger than that 



of the right nanodisk tr( 



— r filter 



T L > Tr. 



(5.12) 



AO = 2\m*L/h, 



(5.10) 



Second, the applied magnetic field is taken so small that the 
nanodisk spin can be controlled by a polarized current. For 
definiteness we take the direction of the magnetic field to be 
up. Third, the lead has the Rashba-type interaction such that 
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(a) without external electric field (spin current flows) 

M M M ♦ M 
I + 

(b) with external electric field (spin current is blocked) 

a) nl CD it * 

FIG. 7: Illustration of spin-field-effect transistor with unpolarized 
current coming from left. We apply the same external magnetic 
field to both of the nanodisks. The unpolarized current is filtered 
by the left nanodisk, and only up- spin electrons go through the cen- 
teral lead, (a) Without external electric field, the electron spin in the 
central lead does not rotate, and the outgoing spin-polarized current 
exists, (b) With an appropriate external electric field, since the elec- 
tron spin in the central lead rotates by the Rashba-type interaction, 
the outgoing current does not exist. 



the rotation angle is AO = tt — S with small 5, say, S ~ 
0.17T. When no currents enter the nanodisk, the direction of 
two nanodisk spins is identical due to a tiny external magnetic 
field. When we inject the current in this state, the net outgoing 
current is very small, 



cos ■ 



0. 



(5.13) 



This is the "off" state of spin diode. 

spin diode 




FIG. 8: Illustration of spin diode made of two nanodisks with dif- 
ferent size. By controlling the bias voltage A/i, the current flows 
from the left lead to the right lead (A/i > 0), or in the opposite way 
(A/i < 0). The incoming current is unpolarized, which is made po- 
larized by the first nanodisk. The electron spin in the central lead is 
rotated by the Rashba-type interaction. 

Let us inject an unpolarized pulse square current to the sys- 
tem, starting at t = t{ and finishing at t = tf , 



i a (t) = r n o(t-t i )o(t f -t), 



(5.14) 



where a denotes the spin. The system become the "on" state 
by the pulse. When the bias voltage is such that A/i > 0, 
the current flows into the left nanodisk and then into the right 
nanodisk. The left nanodisk acts as a spin filter. The current 
in the central lead is initially up-polarized but is rotated by 



the angle AO due to the Rashba-type coupling effect. Then it 
enters the right nanodisk. Since the relaxation time is r R , the 
total spin-dependent charge carried by the current is given by 



Q T = / J t <ft = Q T (t f ) - Q T (U) , 



Qi= J * ' hdt = Q i (tf)-Qi (U) 



with 



Q^ (t) =I in r R sinrT 1 [exp [(t - t ) /r R ]] , 
Q i (t) =P n r R shirr 1 [exp [- (t - t ) /r R ]] . 



(5.15a) 
(5.15b) 



(5.16a) 
(5.16b) 



On the other hand, when A/i < 0, the current enter the right 
nanodisk and goes out from the left nanodisk. Since the relax- 
ation time is tl, the total spin-dependent charge is given by 
the above formulas but with the replacement of tr by tl. Be- 
cause the size of two nanodisk are different, these two currents 
behave in a different way. 



1.0 
c 0.8 
fc 0.6 
3 0.4 

0.2 



(a) incoming pulse 



1.0 
0.8 
0.6 
0.4 
0.2 



5 10 15 

time (t/T) 



20 



(b) outgoing wave 



L-»R f 









5 10 15 

time (t/T) 



20 



FIG. 9: Pulse wave. The horizontal axis is the time and the vertical 
axis is the current, (a) The incoming pulse wave, (b) The outgoing 
waves for L — > R and R — > L. They are very different for the same 
incoming pulse wave due to the difference in the nanodisk size. We 
have set U = 5, tf = 15, and tl = 3tr. 

We define the direction dependent total charge 
Q a (L -> R) and Q a (R -> L), where Q a (L -> R) is 
the total charge with the spin a when charges flow from left 
to right, while Q a (R — > L) is the total charge when charges 
flow from right to left. We find the relation 

Q T (L^R)> Q t (R -> L) » Q i (R^L)> Q i (L - R) 

(5.17) 

from (15.12b , which implies the up (down) component in- 
creases (decreases) from the initial value. 
We define the rectification coefficient by 



R 



Qr (L -+ R) + Qi (L - R) 



which is approximately equals to 



L) 



R: 



(5.18) 



(5.19) 



We illustrate this rectification coefficient as a function of the 
pulse width At = tf — U for various initial phases 9(0) or 
equivalently to in Figll0| Each curve has a peak structure. 
Hence, when the relaxation-time ratio tl/tr and the initial 
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FIG. 10: The vertical axis is the rectification coefficient R. The hor- 
izontal axis is the pulse width At = tf — U. We plot R for various 
initial values £o/tr = log tan [(tt — 5) /2]. We have set tl = 3tr. 

phase (0) are given, it is possible to optimize the width At 
so that the rectification coefficient R is maximized. This max- 
imized value of R diverges as 6 (0) — > 7r. However, the relax- 
ation time diverges as well. It would be efficient to take the 
initial phase (0) to make £ /tr — 5 for a spin diode. 



VI. CONCLUSIONS 



We have studied the electromagnetic properties of the 
zigzag trigonal nanodisk by projecting the system to the zero- 
energy sector. We may regard it as a quasi-ferromagnet char- 
acterized by the exchange energy as large as the Coulomb 
energy. The system is well approximated by the SU(N) 
Heisenberg-Hubbard model. The relaxation time is finite but 
quite large even if the size is very small. Being a ferromagnet, 
it can be used as a spin filter. Namely, only electrons with spin 
parallel to the spin of the nanodisk can go through it. Addi- 
tionally, it has a novel feature that it is not a rigid ferromagnet. 
The incoming spin-polarized current can rotate the nanodisk 
spin itself. Combining the advantages of both these properties, 
we have proposed a rich variety of spintronic devices, such as 
spin memory, spin amplifier, spin valve, spin-field-effect tran- 
sistor and spin diode. Graphene nanodisks could well be basic 
components of future nanoelectronic and spintronic devices. 

I am very much grateful to N. Nagaosa and S. Tarucha for 
many fruitful discussions on the subject. 
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